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Abstract 

Relativistic Hartree equations for spherical nuclei have been derived from a rel- 
ativists quark model of the structure of bound nucleons which interact through the 
(self-consistent) exchange of scalar (a) and vector (u and p) mesons. The coupling 
constants and the mass of the a-meson are determined from the properties of sym- 
metric nuclear matter and the rms charge radius in 40 Ca. Calculated properties of 
static, closed-shell nuclei from 16 O to 208 Pb are compared with experimental data 
and with results of Quantum Hadrodynamics (QHD). The dependence of the results 
on the nucleon size and the quark mass is investigated. Several possible extensions 
of the model are also discussed. 
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1 Introduction 



Do quarks play an important role in finite nuclei ? This is one of the central questions 
in nuclear physics. We now know that explicit quark degrees of freedom are certainly 
needed to understand deep-inelastic scattering at momentum transfers of several GeV, 
e.g., the nuclear EMC effect [|]. However, it is perhaps not surprising that such physics 
seems irrelevant at a scale of a few tens of MeV, and the idea that the quark structure 
of the nucleon should be of any importance in nuclear physics is often dismissed out of 
hand. 

Within the traditional approach to nuclear physics these tens of MeV are the result of 
a fairly fine-tuned cancellation between short-range repulsion, fairly strong intermediate 
range attraction and large kinetic energy. Nearest neighbour nucleons in nuclear matter 
tend to be an average 1.8 fm apart, which corresponds to the intermediate range attraction 
associated with two-pion-exchange. This attraction is mainly generated by the coupling 
to the A (through the process NN — > NA — > NN), which is just the first excited state 
of the internal quark structure of the nucleon, namely spin-flip. Of course, one can treat 
the NN — > NA system in a coupled channels formalism without ever mentioning quark 
structure. One can even replace the NA coupling by simple scalar meson (a) exchange 
without significantly altering the binding energy of nuclear systems. If this were the full 
story quarks would be irrelevant. 

However, it is well known that two-body interactions alone cannot reproduce either 
the saturation properties of nuclear matter or the binding energies and charge radii of 
finite nuclei. One needs some form of three-body force, and this leads to a nightmare of 
badly known form-factors and coupling constants, of inconsistent treatments of the 7rN, 
NN and NNN systems, and so on||. 

Conventional descriptions of finite nuclei are usually based on the Schrodinger equation 
which involves nucleons interacting through static potentials. These effective potentials 
are sometimes determined partly by the Brueckner-Hartree-Fock equations, and often in- 
volve some phenomenological dependence on the nuclear density to include the effect of 
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higher-order contributions to the N-N force in nuclei [H]. Those density-dependent poten- 
tials are however purely phenomenological. Although those conventional approaches give 
rather successful results for the ground states of nuclei, discrepancies between experimen- 
tal data and theoretical calculations still remain. 

Furthermore, there are several problems at a low energy scale, which might not be ex- 
plained by the traditional approach. One example is the Okamoto-Nolen-Schiffer (ONS) 
anomaly @], which is a well known, long-standing problem in nuclear physics. Conven- 
tional nuclear contributions to the anomaly are thought to be at the few percent level and 
cannot explain the experimental findings^]. The effects of charge symmetry breaking in 
the nuclear force, especially p-u mixing, seem to reduce the discrepancy [|], |5|, 0. However, 
recent investigations of the off-shell variation of the p-u mixing amplitude have put this 
explanation into question 0. 

Another example concerns the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements. There is also a well known discrepancy between the value of the Fermi 
decay constant extracted from super-allowed Fermi beta-decay of nuclear isotriplets and 
that required by unitarity of the CKM matrix. This discrepancy remains at the level of a 
few tenths of a percent (but beyond two standard deviations from the experimental data) 
after the most rigorous investigation of conventional nuclear and radiative corrections ||. 

These facts seem to indicate that the traditional approach may have its limitations 
and suggest a need for the study of alternate approaches including sub-nucleonic degrees 
of freedom. 

Exciting a single quark in the nucleon costs about 400 MeV. This is not significantly 
different from the energy required to excite a A. It is the same order of magnitude as 
the scalar and vector potentials required in Quantum Hadrodynamics (QHD)[||. Further- 
more, the quarks are very light and should be able to respond rapidly to their environment. 
We know of no physical argument why this response should be ignored. 

About eight years ago Guichon[|10[ proposed an entirely different model for nuclear 
matter, based on a mean-field description of nucleons described by the non-overlapping 
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MIT bag model, namely the quark-meson coupling (QMC) model. Later this model was 



refined by Saito and Thomasjll] who clarified its relationship to QHD. Related work has 



been carried out by a number of groups [jjj ] . 

The QMC model may be viewed as an extension of QHD in which the nucleons still 
interact through the exchange of a and u mesons. However, the mesons couple not to 
point-like nucleons but to confined quarks. In studies of infinite nuclear matter [jlO, |ll 



it was found that the extra degrees of freedom provided by the internal structure of the 
nucleon mean that one gets quite an acceptable value for the incompressibility once the 
coupling constants are chosen to reproduce the correct saturation energy and density 
for symmetric nuclear matter. This is a significant improvement on QHD at the same 
level of sophistication. A wonderful feature of this picture is that most specific details 
of the model are irrelevant. For example, it does not matter that one uses the MIT bag 
model, all that is needed is that the quarks are relativistic. It does not matter that the u> 
generates the repulsion, all that is needed is that it has a vector character. 

There have been several interesting applications to the properties of finite nuclei using 
the local-density approximation. Surprisingly, the model can provide a semi- quantitative 



explanation of the ONS anomaly when quark mass differences are included fl3|, H| . An 
application of the model, including quark mass differences, has also suggested a previously 
unknown correction to the extraction of the CKM matrix element, V U d, from super-allowed 
Fermi beta-decay, which would bring the discrepancy in the unitarity problem of the CKM 
matrix back to only one standard deviation fl5|, |16|| . Furthermore, in the light of current 
experimental work in relativistic heavy ion collisions [[17], [Tj|, which produce nuclear matter 
at densities several times normal, there has been some initial work on the variation of 
baryon and meson properties with density Finally the model has been applied to 
the case where quark degrees of freedom are undisputedly involved - namely the nuclear 
EMC effect 0. 

However, the inherent problems of the local-density approximation mean that these 
applications can at best be semi- quantitative and it was clearly very important that the 
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extension to finite nuclei be developed. Recently we have succeeded in developing a 



formulation of the quark- meson coupling model for finite nuclei [21], based on the Born- 
Oppenheimer approximation. The effective equation of motion for the nucleon, as well 
as self-consistent equations for the meson mean fields have been derived. We have also 
shown some initial results for 16 and 40 Ca, and discussed, in particular, the spin-orbit 
force in the model and its relation to the corresponding force in conventional models 
involving meson exchange between point-like nucleons. Related work has been carried 
out by several groups [^2], 

Our aim in this paper is to show details of our calculations for the properties of 
spherically symmetric, closed- shell nuclei from 16 O to 208 Pb. 

In Sect .El we briefly summarize the QMC model for both infinite nuclear matter and 



finite nuclei. Starting from the Lagrangian density proposed in Ref. |2jJ we derive the 
Dirac equation for the nucleon and the equations for the meson fields in mean-field ap- 
proximation and discuss the effective cr-N coupling constant. The properties of infinite 
nuclear matter are then reviewed briefly in Sect jO] . We introduce the couplings of the 
p meson and the photon to the quarks, in Sect.2 - 2, in order to describe realistic finite 
nuclei, and again derive explicit forms for the equations of motion of the nucleon and the 
mesons. In Sect.|3], the parameters in our calculations are determined so as to reproduce 
the rms charge radius of 40 Ca, and then our calculated results are compared with those 
of QHD and the experimental data. Our results agree with data favorably. We give a 
summary and discuss some remaining problems and some possible extensions in SectM 



2 Relativistic formulation of the QMC model 

The solution of the general problem of a composite, quantum particle moving in back- 
ground scalar and vector fields that vary with position is extremely difficult. One has a 
chance to solve the particular problem of interest to us, namely light quarks confined in 
a nucleon which is itself bound in a finite nucleus, only because the nucleon motion is 
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relatively slow and the quarks highly relativistic. Thus the Born-Oppenheimer approxi- 
mation, in which the nucleon internal structure has time to adjust to the local fields, is 
naturally suited to the problem. It is relatively easy to establish that the method should 



be reliable at the level of a few percent |2T ] . 

Even within the Born-Oppenheimer approximation, the nuclear surface gives rise to 
external fields that may vary appreciably across the finite size of the nucleon. In Ref . |2l| , 
our approach was to start with a classical nucleon and to allow its internal structure (quark 
wavefunctions and bag radius) to adjust to minimise the energy of three quarks in the 
ground-state of a system consisting of the bag plus constant scalar and vector fields, with 
values equal to those at the centre of the nucleon. Of course, the major problem with the 
MIT bag (as with many other relativistic models of nucleon structure) is that it is difficult 
to boost. We therefore solve the bag equations in the instantaneous rest frame (IRF) of 
the nucleon - using a standard Lorentz transformation to find the energy and momentum 
of the classical nucleon bag in the nuclear rest frame. Having solved the problem using 
the meson fields at the centre of the "nucleon" (which is a quasi-particle with nucleon 
quantum numbers), one can use perturbation theory to correct for the variation of the 
scalar and vector fields across the nucleon bag. In first order perturbation theory only the 
spatial components of the vector potential give a non- vanishing contribution. (Note that, 
although in the nuclear rest frame only the time component of the vector field is non-zero, 
in the IRF of the nucleon there are also non- vanishing spatial components.) This extra 
term is a correction to the spin-orbit force, which is however very small[^TJ. Blunden 



and Miller [23[ have also investigated the QMC model for finite nuclei recently. They too 
reported that the correction due to the variation of meson fields inside the nucleon is not 
large. 



As shown in Ref. [f2lj , the basic result in the QMC model is that, in the scalar (a) and 
vector (u) meson fields, the nucleon behaves essentially as a point-like particle with an 
effective mass Mjy-, which depends on the position through only the a field, moving in a 
vector potential generated by the u meson. Because of the vector character, the vector 
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interactions have no effect on the nucleon structure except for an overall phase in the wave 
function, which gives a shift in the nucleon energy. 

If we restrict ourselves to consider static, spherically symmetric nuclei, the effective 
Lagrangian density involving the quark degrees of freedom and the meson fields in mean- 
field approximation would be given by [f2l[ 



£qmc = #7 • d - M^(a(r)) - g^f)^ 

~ ^[(Va(rl) 2 + mM^ + ^[(Vc,(rl) 2 + m>(f) 2 ], (1) 

where ip(f), &(f) an d 00(f) are respectively the nucleon, a and the time component of 
10 fields. (In the mean field approximation, the meson fields are given by their time 
independent expectation values in the ground state of the nucleus. For a symmetric, 
parity-eigenstate, the spatial component of the 00 field vanishes.) m a and m w are respec- 
tively the masses of the a and 00 mesons. g w is the u-N coupling constant (which is related 
to the corresponding quark-cj coupling constant, g^, as g u = 3g,?,fll0|, [11] , pT|). 



The effective nucleon mass is given by a model describing the nucleon structure 
and it depends on the position through only the a field (as mentioned above). We use 
the MIT bag model in this paper. (A relativistic oscillator model has also been used in 
Ref.pB"] as an alternative model.) If we define the field-dependent cr-N coupling constant, 

M* N {*(r)) =M N - g a (a(r))a(f), (2) 

where M N is the free nucleon mass, it is easy to compare with QHD||. The Lagrangian 
density then becomes 

£qmc = ■ d - M N + g a (a(r))o-(r) - g w uj(r)-f ]ip 

- ^[(Va(rl) 2 + m>(rl 2 ] + i[(Vc,(rl) 2 + m>(r) 2 ]. (3) 

The difference from QHD clearly lies in the fact that the internal structure of the nucleon 
has forced a known dependence of the <r-N coupling constant on the scalar field itself. 
(Note that this dependence is not the same as that adopted in the density-dependent 
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hadron field theory p4f where the meson-nucleon vertices are assumed to depend directly 
on the baryonic densities.) 

Equation ([!]) yields the Dirac equation for the nucleon 

[i 7 • d- M^{a)- ^70^ = 0, (4) 

as well as the equations for the meson mean-fields 

(-V r 2 +m>(r) = -(J^M* N (<7))ps(7), (5) 
(-V 2 + ml)u(r) = ^p B (f), (6) 

where p s and p B are respectively the scalar and baryon densities of the nucleon in the 
nucleus, A, which are defined by the expectation values in the ground state of A: 

Ps (r) = (A\^(f)\A) and p B (f) = (A\tfi/;(f)\A). (7) 

One can easily see that the derivative of Mjy with respect to a in Eq.(|J) is the response 
of the nucleon to the external scalar field, and that it is given by the scalar density of a 



quark in the nucleon [[11 



'dM* N 



-Zglf driP q iP q = -3glS(a), (8) 

J ban 



\ da J ° Jbag 

where g% is the quark-a coupling constant and ip q is the quark wave function in the 
nucleon bag with radius R* B in matter (we denote the bag radius of the free nucleon by 
Rb)- We define S(cr) by Eq.(||), which is a function of scalar field itself. Using the MIT 
bag model|5T[, S is explicitly given by 

S{a{r)) = m-l)+m q R B ^ (9) 



where Vt = ^Jx 2 + (R B m*) 2 is the kinetic energy of the quark and m*(r) is the effective 
quark mass defined by m q — g^cx(f), with the quark mass m q inside the free nucleon. 
Furthermore, we define the scalar density ratio, S(cr)/S(0), by C(cr) and the cr-N coupling 
constant at a = by g a (i.e., g a = g a {a - 



C(a)=S(a)/S(0) and ^ = 3^5(0). (10) 



§ Note the change in notation from the earlier works of Saito and Thomas where C was used for what 
we now call S. 

9 



Comparing with Eq.(H), we find 



, 0n , 9*C(*) = -^\g„(*)*], (11) 
and the equation for the a field becomes 

(-V r 2 + ml)a(r) = g a C{a)p s {r). (12) 
We will discuss the cr-N coupling constant further in the next subsection. 

2.1 Nuclear matter limit 

In isospin-symmetric, infinite nuclear matter the sources of the fields are constant and 
can be related to the nucleon Fermi momentum fcir|J: 

4 r - 2k 3 

<"> = p^i /*»<*' = (13) 

4 r -> M* 
P' = 7^3 dk9(k F -k) N ^ , (14) 



AT 

where Mjy means the constant value of the effective nucleon mass given by the MIT bag 
model (detailed derivations are given in Ref. ||21|). 

As in Ref.||21||, the bag constant B and the parameter zq (which accounts for the sum 
of the cm. and gluon fluctuation corrections) in the familiar form of the MIT bag model 
Lagrangian are fixed to reproduce the free nucleon mass (M^ = 939 MeV) under the 
stationary condition in terms of the free bag radius, Rb- In the following we treat the 
free bag radius as a variable parameter to test the sensitivity of our results to the free 
nucleon size. The results for B and Zq are shown in Table p]. We choose the free quark 
mass m q = 0, 5, 10 MeV. 

Let (a, uj) be the constant mean- values of the meson fields. From Eqs.(^) and (12) we 
find 

v = — 5~> ( 15 ) 



m 2 



* 



9*„,^ 4 f jr.nn. ^ M N 



W = S^T^s dk9(k F -k) « , (16) 
m 2 a (2tt)3 J / M *2 + k2 
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Table V. Bag constant and parameter z for several bag radii and quark masses. 



m 9 (MeV) 





5 


10 




0.6 0.8 1.0 


0.6 0.8 1.0 


0.6 0.8 1.0 


B 1 /4(MeV) 
z 


211.3 170.3 144.1 
3.987 3.273 2.559 


210.9 170.0 143.8 
4.003 3.295 2.587 


210.5 169.6 143.5 
4.020 3.317 2.614 



where C(a) is now the constant value of C in the scalar field. This self-consistency 
equation for a is the same as that in QHD, except that in the latter model one has 



C(a) = 1 [|TT| p3| . (This corresponds to the heavy-quark mass limit of the QMC model. 

950 i ^ 1 1 — 



> 




2.5 



Pb'Po 



Figure 1: Effective nucleon mass as a function of density (for m g =5MeV). The solid, 
dotted and dashed curves correspond to Rb = 0.6, 0.8 and 1.0 fm, respectively. 

Once the self-consistency equation for a has been solved, one can evaluate the total 
energy per nucleon [|K], [U], f2~T|j : 

4 



j^jtotal I 



mla 2 



glpB 



, , dkd{k F - k)TjM# + k 2 + -2— + , 

[2'kYpb J 2p B 2ml 



(17) 



We determine the coupling constants, g a and g u , so as to fit the binding energy (—15.7 
MeV) per nucleon and the saturation density, po — 0.15 fm -3 {k% = 1.305 fm -1 ), for 
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symmetric nuclear matter. The coupling constants and some calculated properties of 
nuclear matter (for m a = 550 MeV and m w = 783 MeV) at the saturation density are 
listed in Table @. The last three columns show the relative changes (from their values 
at zero density) of the bag radius (5R* b /Rb), the lowest eigenvalue (5x/x ) and the rms 
radius of the nucleon calculated by the quark wave function (Sr*/r q ) at saturation density 

Table 2: Coupling constants and calculated properties for symmetric nuclear matter at 
normal nuclear density. The effective nucleon mass, Mjy-, and the nuclear incompressibil- 
ity, K, are quoted in MeV. The bottom row is for QHD. 



m g (MeV) 


R B (fm) 


9l/^ 


til** 


M% 


K 


5R B /R B 


5x/x 


5r*/r q 




0.6 


5.84 


6.29 


730 


293 


-0.02 


-0.13 


0.01 





0.8 


5.38 


5.26 


756 


278 


-0.02 


-0.16 


0.02 




1.0 


5.04 


4.50 


774 


266 


-0.02 


-0.19 


0.02 




0.6 


5.86 


6.34 


729 


295 


-0.02 


-0.13 


0.01 


5 


0.8 


5.40 


5.31 


754 


280 


-0.02 


-0.16 


0.02 




1.0 


5.07 


4.56 


773 


267 


-0.02 


-0.19 


0.02 




0.6 


5.87 


6.37 


728 


295 


-0.02 


-0.13 


0.02 


10 


0.8 


5.42 


5.36 


753 


281 


-0.02 


-0.16 


0.02 




1.0 


5.09 


4.62 


772 


269 


-0.02 


-0.18 


0.03 


QHD 




7.29 


10.8 


522 


540 









The most notable fact is that the calculated incompressibility, K, is well within the 
experimental range: K = 200 ~ 300 MeV. Also our effective nucleon mass is much larger 
than in the case of QHD. Although the bag radius shrinks a little at finite density, the 
rms radius of the quark wavefunction actually increases by a few percent at saturation 
density. 
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In Fig.[I], we show the effective nucleon mass in medium, as a function of density, for 
m q = 5 MeV. Its dependence on the quark mass is weak. The scalar density ratio C(a) 



decreases linearly (to a very good approximation) with (yf CT a|^TJ. It is very useful to have 
a simple parametrization for C and the form: 

C(W) = 1 - a x {g a a), (18) 

with g a a in MeV (recall g a = g a {o = 0)) and a = (6.6, 8.8, 11) xlO -4 for R B = (0.6, 
0.8, 1.0) fm, respectively, is quite accurate up to 2 ~ 3p - (These values for the slope 
parameter, a, are valid for all quark masses listed.) 

As a practical matter, it is easy to solve Eq.(|TT|) for g a {(j) in the case where C(a) is 
linear in g a a as Eq.(lll). In fact, we find 



M* N = M N - g a 



1 - \{9o°) 



a. (19) 



The effective <r-N coupling constant, ^^(c"), decreases at half the rate of C(cx). Eq.(|19|) is 
quite accurate up to twice nuclear matter density. 

Having explicitly solved the nuclear matter problem by self-consistently solving for the 
quark wave functions in the bag in the mean scalar field, one can solve for the properties 
of finite nuclei without explicit reference to the internal structure of the nucleon. All one 
needs is Eqs.(|l8D and (pi]) for C(cr) and ilf^ as a function of g a a. 
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2.2 Finite nuclei 

To describe a nucleus with different numbers of protons and neutrons [Z ^ N), it is 
necessary to consider the contributions of the p meson. Any realistic treatment of nuclear 
structure also requires that one introduces the Coulomb force. The interaction Lagrangian 
density for the p- and 7-quark couplings in the nucleon IRF is given by 

£? +7 = -gtfjwYaJRF - eq'{\ + f (20) 

where q is the quark field (the prime means the IRF) and p^ a IRF and A IRF are respectively 
the p meson field with isospin component a and the photon field, is the isospin 

operator for the quarks. g q and e are the quark-p coupling constant and the electric 
charge. 

We saw in Ref. that this leads to both a central and a spin-orbit potential for the 
nucleon in the nuclear rest frame. For the isoscalar, u meson these potentials are well 
represented by a vector tu-nucleon coupling because the isoscalar magnetic moment of the 
nucleon is near unity. However, for the p, the relativistic formulation at the nucleon level 
requires a strong tensor coupling (~ o~^ v q v ) if it is to reproduce the interaction given by 



the quark model[21 



In Hartree approximation a tensor coupling gives only a spin-orbit force for a nucleon 
bound in a static, spherical nucleus. However, in Hartree-Fock it can also give rise to 
a central force which would be necessary to reproduce the bulk symmetry energy. For 
simplicity, in the present work we follow the usual procedure adopted in the Hartree 
treatment of QHD of using only a vector coupling for the p, with g p adjusted to give the 
bulk symmetry energy in MFA - where g p = g q p . (The accuracy of this approximation will 
be examined in detail in future work.) 

In this case, in mean field approximation, only a = 3 contributes in Eq.(^). If we 
denote by b(f) and A(r) the mean values of the time component of p field and the Coulomb 
field in the nuclear rest frame, we can transpose our results for the u field by trivial isospin 
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factors: 



3 ™ ft i ~N\ 



3<?>(0 - 9 P ^-b(r) or -(1 + r 3 w )A(f), (21) 



where is the third component of the nucleon isospin operator. 

In summary, our effective Lagrangian density in mean field approximation takes the 
form: 

£qmc = ^[il ■ d - M N + g a (a(r))a(r) - g u u(r)^ 
- 9 P \Kr) lQ --{l + T?)A{r) 1 ^ 
~ ^(Va(f)) 2 + mla{rf\ + \[{Vuj{r)) 2 + m>(rf] 
+ i[(V6(f)) 2 + m^(r) 2 ] + ^(VA(r)) 2 . (22) 

The variation of the Lagrangian results in the following equations for static, spherically 

symmetric nuclei (see also Eqs.(^) and (0)): 
d 2 2d 

-r^cr(r) + -—<r{r) - m\o{r) = -g a C{a{r))p s {r) 
dr z r dr 

occ 

= -9,C(a(r))Y,d a (r)(\G a (r)\ 2 -\F a (r)\ 2 ), (23) 

a 

d 2 2d 



-9»E d <*(r)(\G a (r)\ 2 + |F Q (r)| 2 ), (24) 



d 2 2 d g 

— b(r) + -—b{r) - m 2 b(r) = -^p 3 {r) 
dr z r dr p 2 



occ 

-fT,Mr)H ta - 1/2 (\G a (r)\ 2 + \F a (r)\ 2 ), (25) 



OCC 1 

-eJ2d a (r)(t a + -)(\G a (r)\ 2 + \F a (r)\ 2 ), (26) 



where d a {r) = (2j a + l)/4irr 2 and 



d k 1 

— G«(r) + -G a (r) - [e Q - g^{r) - t«^6(r) - (t a + -)eA(r) + M N 

dr r 2 

g*Wr))a(r)]F a (r) = 0, (27) 

d k 1 

—F a (r)--F a (r) + [e a -guU}(r)-t a g p b(r) - (t a + -)eA(r) - M N 

+ 9(T (a(r))a(r)}G a (r) = 0. (28) 
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Here iG a (r)/r and —F a (r)/r are respectively the radial part of the upper and lower 
components of the solution to the Dirac equation for the nucleon (a labelling the quantum 
numbers and e a being the energy) under the normalization condition: 

Jdr(\G a (r)\ 2 + \F a (r)\ 2 ) = l. (29) 

As usual, k specifies the angular quantum numbers and t a the eigenvalue of the isospin 
operator C(a) and g a (a) are given by Eqs.@ and (|2|), respectively - or practically, 

by Eqs.(H) and (0), i.e., 



1 - 2^ a \^ 



The total energy of the system is then given by 

OCC « 

E tot = ^(2j a + l)e Q - - J df i-g a C{a(r))a(r)p s (r) 



(30) 



+ 9uu{r)p B {r) + ^g p b(r)p 3 (r) + eA(r)p p (r)}. (31) 



3 Numerical results for finite nuclei 

Equations (^) to ( f29|) give a set of coupled non-linear differential equations, which may be 
solved by a standard iteration procedure [^5j. In this paper the numerical calculation was 
carried out by modifying the technique described by Horowitz et al. [53], . Changing the 



initial condition for the meson fields, the <r-N coupling constant and the scalar density in 
their fortran program[p26[, the calculation can be performed very easily^. The calculation 



is achieved in at most 20 iterations when it is performed with a maximum radius of 12 
(15) fm on a mesh of 0.04 fm for medium mass (Pb) nuclei. 



^If the numerical convergence is slow it may be improved by mixing appropriately the meson potentials 
given by the i-th iteration and those by the (i — l)-th iteration - as is usually done in non-relativistic 
calculations. 
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3.1 Determination of parameters 

There are seven parameters to be determined: g a , g^, g p , e, m a , m u and m p . We take 
the experimental values: = 783 MeV, m p = 770 MeV and e 2 /47r = 1/137.036. The 
coupling constants g a and g u are fixed to describe the nuclear matter properties with m a 
= 550 MeV in Sec.|3 (see Table |). 

The cr-meson mass however determines the range of the attractive interaction and 
changes in m a affect the nuclear-surface slope and its thickness. Therefore, as in the 
paper of Horowitz and Serot[^], we adjust m a to fit the rms charge radius of 40 Ca, 
Tc/i ( 40 Ca) = 3.48 fm, the experimental value |27|. We should notice here that variations of 



m a at fixed {ga/^a) have no effect on the infinite nuclear matter properties. Therefore, 
keeping the ratio (g a /rn a ) constant we vary m a to fit the rms charge radius of 40 Ca. 
We expect that m a ranges around 400 ~ 550 MeV[^3|, The last parameter, g P) is 



adjusted to yield the bulk symmetry energy per baryon of 35 MeV||25||. We summarize 
the parameters in Table H. 



Table 3: Model parameters for finite nuclei. 



m 9 (MeV) 


i2fl(fm) 


9ll^ 


9l/^ 




m a (MeV) 




0.6 


3.55 


6.29 


6.79 


429 





0.8 


2.94 


5.26 


6.93 


407 




1.0 


2.51 


4.50 


7.03 


388 




0.6 


3.68 


6.34 


6.78 


436 
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3.12 


5.31 


6.93 


418 




1.0 


2.69 


4.56 


7.02 


401 




0.6 


3.81 


6.37 


6.78 


443 


10 


0.8 


3.28 


5.36 


6.92 


428 




1.0 
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4.62 


7.02 


416 
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Figure 2: Model predictions for the effective nucleon mass, the baryon and scalar densities 
in 208 Pb (for m q = 5 MeV and Rb = 0.8 fm). The scale on the right vertical axis is for 
M* N /M N . 

In Fig.^|, we show the baryon and scalar densities as well as the effective nucleon mass 
in lead. We expect that the baryon density in the interior of lead would be close to 
the saturation density of infinite nuclear matter. As seen in the figure, the calculated 
baryon density at the center is quite close to 0.15 fm -1 , which supports our choice of the 
parameters. 
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3.2 Comparison with experimental data 

First we show calculated charge density distributions, p c h, in comparison with those of 
QHD||, |25| and the experimental data in Figs.[3]— 112]. Having solved Eqs.(|23|) ~ (^), we 
obtain the point-proton and neutron densities in a nucleus. In addition to those densities, 
we should include the effects of cm. corrections, nucleon form factors, meson-exchange 
currents, etc. However, one knows that the dominant correction comes from the proton- 
form factor p5|. We calculate the charge density as a convolution of the point-proton 
density, p p (r), with the proton charge distribution, p^. h (f): 



P*(f) = J df' p p ch (f- f')p p (f'), (32) 

where we have used a gaussian form for p^, h 

P^(r) = (/5/vr) 3 / 2 exp(-/5r 2 ). (33) 

The parameter /3, which determines the proton size, is chosen so as to reproduce the 
experimental rms charge radius of the proton, 0.82 fm (i.e., (3 = 2.231 fm -2 ). As we noted 
before, in the present model the rms radius of the nucleon in nuclear matter increases 
a little. However, this amount is quite small and it can be ignored in the numerical 
calculations of nuclear parameters. 

The charge density distributions for 40 Ca (for m q = 0, 5, 10 MeV and Rb = 0.8 fm) 
are presented in Figs. |3|— 15|. The experimental data is taken from Ref. ||29||. Once the rms 



charge radius of 40 Ca is fitted by adjusting m a the QMC model can reproduce p c h( 40 Ca) 
quite well. As seen in the figures, the calculated charge densities lie almost within the 
experimental (hatched) area. We note that the dependence of p c h( 40 Ca) on the bag radius 
is quite weak. 

It is interesting to see the quantum oscillations of the interior density in lead. The 
dependence of p c /i( 208 Pb) on the quark mass is illustrated in Figs.|— 1[ We note that 
the dependence of p c / l ( 208 Pb) on the free nucleon size is also quite weak. As seen in the 
figures, our model gives charge densities very close to those of QHD and still somewhat 
larger in the central region than those observed experimentally |3T 
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0.10 




Figure 3: Charge density distribution for 40 Ca (for m q = MeV and Rb = 0.8 fm) 
compared with the experimental data and that of QHD. 



In Figs.y and |H| we show respectively the charge density distributions for 16 O and 
90 Zr. For zirconium the calculated p c h lies in between that of the non-relativistic density- 
dependent Hartree-Fock calculations and that of QHD. The experimental data for 
oxygen and zirconium are taken from Refs. [31]] and [|32|| . (For both cases the dependence 



of p c h on m q and Rb is weak.) 

To see the isotope shift in charge density we have plotted r 2 times the difference 
between p c /,,( 40 Ca) and /^(^Ca) in Figs.|ll| and [12|. Its dependence on the bag radius is 
again weak for a small quark mass while it becomes a little stronger for m q = 10 MeV, as 
shown in Fig.^. The experimental data is taken from Ref.|33[. (Note that in this case 



we also checked that including the charge distribution of the neutron had a small effect.) 

In Figs.|l^— [H| we present the point-neutron density distributions, p n , in calcium and 
lead. For 40 Ca, since the dependence of p n on m q and Rb is again fairly weak, only the 
result for m q = 5 MeV and Rb = 0.8 fm is shown, together with the empirical fit [|34|j 
to proton scattering data. We again find reasonable agreement with the data. For 
the isotope shift of p n ( 48 Ca)— p n ( 40 Ca), the calculated difference is closer to those of non- 
relativistic results than to those of QHD. The neutron density distribution in lead is shown 



20 



0.10 



0.08 



^ 0.06 



qL 0.04 



0.02 - 



0.00 



Walecka—Serot 

experiment 

this work 

(m = 5 MeV, R = 0.8 fm) 




Figure 4: Same as Fig.5 (for m q = 5 MeV and Rb = 0.8 fm). 

in Fig.|TJ| Its behavior is again similar to that of QHD. 

Table [| gives a summary of the calculated binding energy per nucleon (E/A), rms 
charge radii and the difference between nuclear rms radii for neutrons and protons (r n — r p ), 
for several closed-shell nuclei. Since the calculated properties do not depend strongly on 
m q and Rb-, we only list the values for the QMC model with m q = 5 MeV and Rb = 



0.8 fm. References for the experimental values can be found in Ref . ||25|| . While there are 
still some discrepancies between the results and data, the present model provides quite 
reasonable results. In particular, it gives much larger binding energies per nucleon than 
those of QHD while still reproducing the rms charge radii for medium and heavy nuclei 
quite well. 

In Figs.|16| and [T7|, the calculated spectra of 40 Ca and 208 Pb are presented. Because 
of the relatively smaller scalar and vector fields in the present model than in QHD, the 
spin-orbit splittings are smaller. The improvement over QHD in the binding energy per 
nucleon comes at the expense of a reduction in the spin-orbit force. We should note that 
there is a strong correlation between the effective nucleon mass and the spin-orbit force. 
This problem is also discussed in Refs. [23| . 



As a test of the sensitivity of the spin-orbit splitting to features of the model, we con- 
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Figure 5: Same as Fig. 5 (for m q = 10 MeV and Rb = 0.8 fm). 

sider the case of a larger quark mass (noting that the Born-Oppenheimer approximation 
requires that the nucleon motion is relatively slow and the quarks are highly relativistic). 
For example, we have calculated the case m q = 300 MeV (and Rb = 0.8 fm) which is 
a typical constituent quark mass. The calculated spectra for m q = 300 MeV are also 
illustrated in Figs.0 and [17] (QMCH). In this case, the various parameters were <7^/47r = 
5.58, g^/4:TT = 8.51 (to satisfy the saturation condition), g 2 p /^ = 6.45, m CT = 497 MeV 
(to fit the rms charge radius of 40 Ca), K = 334 MeV and M* N = 674 MeV at saturation 
density. We should record that the bag radius in this case increases by 7% at saturation 
density (which is uncomfortably large). The slope parameter in Eq.flTSD is a = 3.9 x 10~ 4 . 
One can expect that a heavy quark mass gives a spectrum closer to those of QHD||rT], |2T| . 
We can see from the figures that the calculated spectra are somewhat closer to the exper- 
imental data. We note that the charge density distributions for 40 Ca and 208 Pb are also 
reproduced well in this case. 
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Figure 6: Charge density distribution for 208 Pb (for m q = MeV and R B = 0.8 fm) 
compared with the experimental data and that of QHD. 




Figure 7: Same as Fig.8 (for m q = 5 MeV and R B = 0.8 fm). 
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Figure 8: Same as Fig. 8 (for m q = 10 MeV and R B = 1.0 fm). 
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Figure 9: Charge density distribution for 16 (for m q = 5 MeV and R B = 0.8 fm) compared 
with the experimental data and that of QHD. 
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Figure 10: Charge density distribution for 90 Zr (for m q = 5 MeV and R B = 0.8 fm) 
compared with the experimental data and that of QHD. 
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Figure 11: Isotope shift between p c /j( 40 Ca) and p c / l ( 48 Ca) compared with the experimental 
data and that of QHD (for m q = and 5 MeV with R B = 0.8 fm). 
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Figure 13: Point-neutron density distribution in 40 Ca (for m q = 5 MeV and R 
fm) compared with that of QHD and the empirical fit. 
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Figure 15: Point-neutron density distribution in 208 Pb (for m q = 5 MeV and Rb = 0. 
fm) compared with that of QHD and the empirical fit. 
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Table 4: Binding energy per nucleon E/A (in MeV), rms charge radius r ch (in fm) and 
difference between r n and r p (in fm). m q = 5 MeV and Rb = 0.8 fm. (* fit) 





—E/A 






Model 


QMC 


QHD 


Exp. 


QMC 


QHD 


Exp. 


QMC 


QHD 


Exp. 


16Q 


5.84 


4.89 


7.98 


2.79 


2.75 


2.73 


-0.03 


-0.03 


0.0 


40 Ca 


7.36 


6.31 


8.45 


3.48* 


3.48* 


3.48 


-0.05 


-0.06 


0.05±0.05 


48 Ca 


7.26 


6.72 


8.57 


3.52 


3.47 


3.47 


0.23 


0.21 


0.2±0.05 


90 Zr 


7.79 


7.02 


8.66 


4.27 


4.26 


4.27 


0.11 


0.10 


0.05±0.1 


208 Pb 


7.25 


6.57 


7.86 


5.49 


5.46 


5.50 


0.26 


0.27 


0.16±0.05 
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Figure 16: Model predictions for the energy spectrum for 40 Ca. QMC(H) denotes 
case for m q = 5 (300) MeV and R B = 0.8 fm. 
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Figure 17: Same as Fig. 18 (for 208 Pb) 



30 



4 Summary, discussion and further applications 

Starting with the quark-meson coupling model, in which quarks confined in nucleon bags 
interact through the exchange of scalar and vector mesons, we have presented a gener- 
alisation of QHD with a a field- dependent scalar coupling to describe finite nuclei quan- 
titatively. The physical origin of this field-dependence, which provides a new saturation 
mechanism for nuclear matter, is the relatively rapid increase of the lower Dirac com- 
ponent of the wavefunction of the confined, light quark. We have then derived a set of 
coupled non-linear differential equations which must be solved self-consistently but which 
are not much more difficult to solve than the relativistic Hartree equations of QHD. Our 
calculations for static, closed-shell nuclei from 16 to 208 Pb reproduce fairly well the 
observed charge density distributions, neutron density distributions etc. 

It will be very interesting to explore the connection between the a field- dependence 
of the variation of the effective cx-N coupling constant, which arises so naturally here, 
and the variation found empirically in earlier work. We note, in particular, that while 
our numerical results depend on the particular model chosen here (namely, the MIT bag 
model), the qualitative features which we find (such as the a field-dependent decrease 
of the scalar coupling etc.) will apply in any model in which the nucleon contains light 
quarks and the attractive N-N force is a Lorentz scalar. Of course, it will be important 
to investigate the degree of variation in the numerical results for other models of nucleon 
structure f23|j. 

In the present model there are, however, still some discrepancies in energy spectra 
of nuclei, in particular, the spin-orbit splittings. In the previous section we have briefly 
discussed one possibility to partly overcome this defect in the present model. As an 



alternative approach, Jin and Jennings |22fl and Blunden and Miller [^J] have considered 
variations of the bag constant B and z parameter in nuclear matter, which have been 
suggested by the fact that quarks are presumably deconfined at high enough densities. 
This was taken to suggest that B might decrease with increasing density. Adopting this 
idea, Blunden and Miller []53| have studied properties of both nuclear matter and finite 
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nuclei. In their approach, B in matter, B*, is given by 



B* = B 



1 - a B - 



(34) 



M N 

where U s is an average scalar potential and a B is an arbitrary parameter. For finite nuclei 
the results move towards those of QHD, and there is an improvement in the spin-orbit 
splittings. However, the nuclear incompressibility, K, is again larger than the experimen- 
tal data, and the bag radius at saturation density increases by 5% ~ 13%[|]| or morc|22|. 



which seems uncomfortably large Q. We probably ought to construct a formalism beyond 
the Hartree approximation for this class of models in the future. 

In this work we have not considered any effect of a decrease of the meson masses 



in nuclear matter. (In Ref.[|19[] these effects were investigated in infinite nuclear matter, 
using an earlier version of the model.) The cx-meson mass which we found it necessary 
to use in this paper ranges from 400 ~ 450 MeV. It seems rather small compared with 
the empirical cx-mass in the OBEP analyses of the free NN force || ||, ||. On the other 



hand, noting its origin in two-pion-exchange with N-A intermediate states, we observe 
that the N-A mass difference would decrease in medium by a similar percentage. It will be 
interesting, in future work, to incorporate a self-consistent reduction of the vector-meson 
mass |T§| (as well as the scalar-meson mass[^B[) into the calculation of the properties of 
finite nuclei. 

The successful generalisation of the QMC model to finite nuclei opens a tremendous 
number of opportunities for further work. For example, earlier results for the Okamoto- 
Nolen-Schiffer anomaly [p^] , the nuclear EMC effect p0[], the charge-symmetry violating 



correction to super-allowed Fermi beta-decay[l5[] and so on, can now be treated in a 



truly quantitative way. Since the present model can provide self-consistent quark wave 
functions in a deeply bound nucleon, e.g., lsi/2 in 208 Pb, it would also be very interesting 
to investigate variations in magnetic moments, Qa and so on, for such nucleons[[37]]. 

For hypernuclei the natural extension (see Ref.JTU]]) is to assume that the a and uj 
mesons couple only to the non-strange constituents. From our discussion of the spin-orbit 



force in Ref.f21|| and the fact that the spin of the A is carried entirely by the strange 
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quark, one can easily see that the A spin-orbit force will arise entirely from the Thomas 
precession term. This means that the A spin-orbit force is very naturally suppressed in 
this model - as observed experimentally. It will be important to follow this observation 
with quantitative results. 

Finally, we list a number of important ways in which this model could be extended: 
replacing the MIT bag by a model respecting PCAC (e.g., the cloudy bag model fl38|l ); 
replacing a-exchange by two-pion exchange; replacing u exchange by nucleon overlap 
at short distance, and so on. In terms of further theoretical development it will also 
be interesting to compare the present model with more phenomenological, non-linear 
extensions of QHD - as reviewed recently in Ref. ||39|| . On the practical side, we stress 
that the present model can be applied to all the problems for which QHD has proven so 
attractive, with very little extra effort. 
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